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Abstract 

A  previous  theory  for  studying  the  distribution  of  non-uniform  fields  in  multiple-quantum-well  photodetectors 
under  an  ac  voltage  is  generalized  by  including  non-adiabatic  space-charge-field  effects.  Numerical  calculations  indicate 
that  field-domain  effects  are  only  important  at  high  temperatures  or  high  voltages  when  both  injection  and  sequential- 
tunneling  currents  are  significant.  On  the  other  hand,  it  is  found  that  the  non-adiabatic  effects  included  in  this  gen¬ 
eralized  theory  become  significant  at  low  temperatures  and  low  voltages  when  field-domain  effects  are  negligible.  In 
order  to  explain  the  non-adiabatic  charge-density  fluctuations  quantum-statistically,  a  non-adiabatic  differential 
equation  is  derived  based  on  the  self-consistent  Hartree  model  by  using  a  shifted  Fermi-Dirac  model  for  the  local 
fluctuation  of  electron  distributions.  The  non-adiabatic  effect  is  found  to  cause  an  “equilibrium”  state  variation  with 
time  under  an  ac  voltage. 
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1.  Introduction 

Multiple-quantum- well  (MQW)  photodetectors 
using  intersubband  transitions  have  attracted  a  lot 
of  studies  over  the  past  few  years  [1].  Transient 
spectroscopy  allows  us  to  gain  information  on  the 
sequential-tunneling  processes  between  quantum 
wells  (QWs)  and  thus  QW  parameters,  including 
geometric  and  QW  capacitances  at  the  same  time. 
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Recently,  a  residual  dark  current  in  quantum-well 
infrared  photodetectors  (QWIPs)  was  reported 
when  an  ac  bias  voltage  was  swept  through  zero 
[2],  Later,  a  roll-off  of  the  dynamical  responsivity 
in  QWIPs  was  observed  when  the  frequency  of  a 
chopped  incident  optical  flux  was  increased  be¬ 
yond  a  certain  value  [3],  More  recently,  a  counter¬ 
clockwise  hysteresis  loop  for  the  tunneling  current 
and  a  clockwise  hysteresis  loop  for  the  emission 
current  in  QWIPs  was  seen  experimentally  as  the 
device  temperature  was  swept  up  from  10  to  300  K 
and  then  back  down  [4],  These  new  phenomena 
found  in  QWIPs  were  physically  explained  by  a 
non-adiabatic  sequential-tunneling  model  [5].  The 
non-adiabatic  effect  discussed  in  this  paper  refers 
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to  the  fact  that  a  transient  current  not  only 
depends  on  the  magnitude  of  an  electric  field,  but 
also  depends  on  its  time  derivative. 

For  the  field-domain  effect,  we  generalize  the 
previous  theory  [6]  for  studying  the  distribution 
of  non-uniform  fields  in  MQW  photodetectors 
under  an  ac  voltage  by  including  non-adiabatic 
space-charge-field  effects.  We  find  from  numerical 
results  that  field-domain  effects  are  only  impor¬ 
tant  at  high  temperatures  or  high  voltages  when 
both  injection  and  sequential-tunneling  currents 
are  significant.  On  the  other  hand,  we  find  that 
non-adiabatic  effects  included  in  the  generalized 
theory  become  visible  at  low  temperatures  and 
low  voltages  when  field-domain  effects  are  negli¬ 
gible.  The  time  duration  for  non-adiabatic  effects 
is  found  to  depend  on  the  quantum  capacitance, 
while  the  classical  dielectric  displacement  current 
is  found  to  be  related  to  the  geometric  capaci¬ 
tance.  Moreover,  a  negative  conduction  current 
is  predicted  under  a  positive  voltage  in  non¬ 
steady  state. 

For  non-adiabatic  transport,  we  consider  elec¬ 
trons  in  a  MQW  system  in  the  presence  of  a  uni¬ 
form  ac  electric  field.  We  assume  that  electrons 
during  the  sequential-tunneling  process  only  see  an 
instantaneous  electric  field  and  stay  in  equilibrium 
states  due  to  very  fast  elastic  and  inelastic  scatter¬ 
ing  inside  the  quantum  well.  The  non-adiabatic 
effect  causes  an  “equilibrium”  state  to  vary  with 
time.  As  a  result,  the  charge-density  fluctuation  in 
the  QW  will  modify  the  Hartree  potential  in  the 
surrounding  barrier  region,  and  thus  greatly  affect 
the  sequential  tunneling  of  electrons  through  the 
barrier.  Simultaneously,  the  non-adiabatic  fluctu¬ 
ation  of  charge  density  also  modifies  electronic 
states  in  the  quantum  well  within  the  self-consistent 
Flartree  model  [7]. 

The  organization  of  the  paper  is  as  follows.  In 
Section  2,  we  introduce  our  model  for  classical 
charge-density  fluctuations  and  the  distribution  of 
field  domains  in  MQWs.  In  Section  3,  we  gener¬ 
alize  the  model  to  include  quantum  non-adiabatic 
charge-density  fluctuations  and  their  effects  on  the 
distribution  of  field  domains  in  MQWs.  Section  4 
is  devoted  to  the  derivation  of  a  non-adiabatic 
differential  equation  based  on  the  self-consistent 
Hartree  model  by  using  a  shifted  Fermi-Dirac 


model  for  the  local  fluctuation  of  electron  dis¬ 
tributions  and  provides  a  quantum-statistical 
explanation  for  the  non-adiabatic  charge-density 
fluctuations  introduced  in  Section  2.  Some 
numerical  results  for  classical  and  quantum 
charge  fluctuations,  as  well  as  for  non-adiabatic 
fluctuations  in  electron  distributions,  are  pre¬ 
sented  in  Section  5.  The  paper  is  concluded  in 
Section  6. 


2.  Classical  charge  fluctuations 

In  this  section,  we  study  effects  of  field  domains 
resulting  from  the  imbalance  between  injection 
and  sequential-tunneling  currents  and  show  that 
these  effects  become  negligible  at  low  temperatures 
and  voltages. 

Electrons  in  QWs  are  confined  in  the  direction 
perpendicular  to  the  wells,  while  electrons  in 
heavily-doped  contact  layers  are  free  in  all  three 
directions.  Therefore,  we  expect  the  tunneling  of 
electrons  from  a  contact  layer  to  a  QW  (3D-to-2D) 
will  be  physically  different  than  that  from  one  QW 
to  another  (2D-to-2D). 

As  shown  in  Fig.  1,  we  see  that  the  distribution 
of  uniform  dc  electric  fields  in  (a)  with 
<?o  =  =  ■  ■  •  =  $n  is  not  stable  if  the  in¬ 

jection  current  flowing  from  the  left  contact  layer 
to  the  first  QW  and  the  sequential-tunneling  cur¬ 
rent  flowing  from  the  first  QW  to  the  second  QW 
are  different  [6,8,9].  As  an  example,  we  assume  in 
(b)  that  the  injection  current  is  smaller  than  the 
sequential-tunneling  current.  In  this  situation,  the 
local  field  S0  in  the  first  (emitter)  barrier  layer  has 
to  be  increased  so  as  to  equalize  these  two  cur¬ 
rents.  As  a  result  of  S0  >  <$\  =  $2  =  •  •  •  =  <%n,  we 
know  from  the  Maxwell  equations  that  the  charge 
density  in  the  first  QW  will  be  reduced  relative  to 
the  others  [8],  i.e.  p{  <  p2  =  p2  =  ••  •  =  pN.  Now, 
let  us  further  compare  the  tunneling  currents 
flowing  from  the  first  QW  to  the  second  QW  and 
that  flowing  from  the  second  QW  to  the  third  QW. 
We  realize  from  (c)  that  the  tunneling  current 
flowing  from  the  first  QW  is  less  than  that  flowing 
from  the  second  QW  since  px  <  p2  for  S\  =  S2. 
Therefore,  S\  >  S‘2  =  ■  ■  ■  =  S’N  is  required  so  as  to 
equalize  these  two  sequential-tunneling  currents. 
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Fig.  1.  The  processes  (a)-(d)  for  the  formation  of  field  domains  in  a  multiple-quantum-well  (MQW)  sample  with  N  quantum  wells 
(QWs).  Here,  Sk  for  k  =  0,1,2, ...  ,N  indicates  local  electric  fields  in  (N  +  1)  different  barrier  layers,  and  pk  for  k  =  1,2, ...  ,N  cor¬ 
responds  to  different  charge  densities  inside  N  QWs.  /(„  is  the  chemical  potential  in  contact  layers,  Ft,  is  the  applied  voltage,  and  L,  is  the 
total  length  of  MQWs. 


Consequently,  we  are  left  with  p{  <  p2  < 
p2  =  '  "  =  Pn-  This  process  will  continue  up  to  the 
last  (ATh)  QW,  as  displayed  in  (d),  until  the  initial 
field  distribution  So  >  S\  >  S2>  ■  ■  ■  >  SN  and 
density  distribution  px  <  p2  <  p3  <  ■  ■  ■  <  pN  are 
reached.  Here,  the  local  field  Sj  is  a  constant  in  the 
/'th  barrier  layer,  and  the  splitting  of  different  local 
electric  fields  constructs  the  field  domains  in 
MQWs.  The  electric  fields  in  different  barrier  re¬ 
gions  in  steady  state  will  be  redistributed  if  the 
total  tunneling  current  flowing  into  the  bottom 
contact  layer  is  different  from  the  injection  current 
flowing  out  of  the  top  contact  layer.  Therefore,  the 
above  field-adjustment  process  will  be  repeated 
again  and  again  until  the  total  tunneling  current 
flowing  to  the  bottom  contact  layer  is  equal  to  the 
injection  current  flowing  from  the  top  contact 
layer,  and  then,  a  stable  distribution  of  local  fields 
is  formed.  It  should  be  noted  that  the  distribution 
of  electric-field  domains  in  steady  state  depends  on 
the  value  of  voltage  applied  at  each  moment.  A 
variation  of  the  applied  voltage  with  time  results  in 
charge-density  oscillations  in  each  QW,  and 
varying  field  domains  as  well. 

The  injection  current  density  (3D-to-2D)  from 
the  top  contact  layer  to  the  first  QW  is  calculated 
to  be  [8] 


WnflhoThWWoW] 


em*kBT  So(t) 


'  2n2k3  |  AW  I 

r 

x  In 


&Efr[E,\So{i)\] 


_ 1  +  exp  [Wo  —  E)/kBT] _ 

1  +  exp  [(/q (t)  ~E-eLB\S0(t)\)/kBT] 


(1) 


where  t  is  the  time,  m*  is  the  effective  mass  of 
electrons,  T  is  the  electron  (or  lattice)  temperature, 
Lb  is  the  thickness  of  the  barrier,  3T\E,  \S(t)\]  is  the 
quantum  transmission  of  electrons  [10]  with  ki¬ 
netic  energy  E  through  a  barrier  biased  by  an 
electric  field  S(t),  and  the  chemical  potential  p0  in 
the  contact  layer  is  related  to  the  electron  con¬ 
centration  nc  by 


«c 


1 

2n2 


3/2 


1  +  exp 


f  E-Po 
V  kB  T 


-i 


For  low  T,  electrons  in  QWs  can  only  populate  the 
ground  subband  with  quantized  energy  E0.  Fur¬ 
thermore,  the  sequential-tunneling  current  density 
(2D-to-2D)  from  the  kth  QW  to  the  neighboring 
(A-+  l)th  QW  is  found  to  be  [1] 
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pm*  r+°° 

=  -I-vdk[Mt)}  /  dE^[E  +  EMm 

Tin  jL\v  Jo 

X  {fo[E-nk(t)\ -f0[E-pk+l{t)+eLB\gk{t)\}}, 

(2) 


where  Tw  is  the  well  width,  A  =  1,2, . . . ,  N  repre¬ 
sents  the  index  of  N  wells,  k  =  N  +  1  corresponds 
to  the  bottom  contact  layer,  pN+l  =  /<0,  f0(x)  = 
[1  +  exp(x/AB7’)]_1,  and  the  chemical  potential 
pk(t)  in  the  Ath  QW  introduced  in  Eqs.  (1)  and  (2) 
can  be  determined  by  the  electron  density  nk(t)  in 
the  Ath  QW  through 


Hk(t)  =  AB7Tn 


exp 


n  h2nk(t) 
m*k\{T 


1 


which  is  measured  from  E0.  In  Eq.  (2),  vk[gk{t)\  is 
the  drift  velocity  of  electrons  in  the  Ath  barrier 
layer,  given  by  [1] 


with  the  saturation  velocity  vs  and  saturation  field 
gs,  respectively. 

As  explained  in  Fig.  1,  the  distribution  of  field 
domains  changes  with  the  value  of  the  voltage 
Vb(t)  applied  to  the  sample.  When  Vb(t)  varies  with 
time,  the  charge  density  will  fluctuate  in  each  QW, 
accompanied  by  field  domains  that  vary  adiabati- 
cally  across  the  whole  sample.  In  this  case,  the 
charge-density  fluctuation  Pj(t')  =  enj(t)  can  be 
described  by  the  equation  [6] 


dPj{t ) 

ck 

E  ( 1  -  rkXjp  [nk{t),nk+1  (0 ,  sk{i)]  e[gk{t)\ 

k=  1 

+ pe(  1  -  PcY-'jg  [ho,  hi  (0 ,  <?»  (01 0[*>(O1 

-Pc  E  (1 

k=j+ 1 

X  0[-^_i(O]^Pc(l -Pc)N~J^\Po,PN(t)M)} 

x  B[-SN(t)]  -Jf  [pj{t)^j+i(t),gj(t)\  0[gj{t)\ 

+  JjD  [h ./ ( ■ 0 , hy- 1  ( t) ,  gj- 1  (i t ) \  6 [  - . gj- 1  (i t) \ , 

(3) 


where  the  small  diffusion  current  [6,9]  is  neglected, 
0(x)  equals  one  for  x  >  0  and  zero  for  x  <  0, 
j  =  1 , 2, . . . ,  N,  and  0  <  Pc  ^  1  is  the  capture 
probability  of  electrons  into  the  QW.  The  first  two 
terms  in  Eq.  (3)  represent  the  forward  contri¬ 
butions  from  capture  current  into  the  /th  QW, 
while  the  third  and  fourth  terms  represent  the 
backward  contributions  from  capture  current  into 
that  QW.  The  last  two  terms  correspond  to  the 
forward  and  backward  tunneling  currents  flowing 
out  of  the  /th  QW,  respectively. 

If  electrons  in  QWs  can  be  simply  viewed  as  a 
distribution  of  sheet  charges  (zero-thickness),  we 
get  the  following  boundary  conditions  [6]  from  the 
Maxwell  equations  for  two  local  fields  on  both 
sides  the  yth  QW 

S j{t)  —  S j-\{t)  — - —  en2n\,  (4) 

£o£r 

where  quantum  properties  of  an  electron  gas  in  a 
QW  have  been  ignored,  m2d  is  the  electron  number 
density  in  equilibrium,  j  =  1, 2, . . . ,  N  and  er  is  the 
relative  dielectric  constant  of  the  well  material. 
Finally,  the  sum  of  individual  voltage  drops  on 
each  period  (well  plus  barrier)  is  fixed  by  the 
voltage  Vb(t).  This  restraint  gives  rise  to 

N 

LagoO)  +  C^b  +  -Aw)  E/  =  Pb{t)  (5) 

k=  1 

and  Lt=LB+  N(LB  +  Tw)  is  the  total  length  of  the 
MQW  structure  in  Fig.  1. 

Combining  Eqs.  (3)  and  (4)  for  forward  con¬ 
tributions  we  know  that  the  sum  of  the  displace¬ 
ment  and  conduction  currents  is  a  constant  for 
Pc  =  1,  i-e. 

£o£r  -  +  jj°i  [fij- 1  {t) ,  Hj(t) ,  gj- 1  (f)] 

=  £0£r  +  JjD  [PjO)  ,  Pj+ 1  (0 ,  ■ 

However,  the  conduction  current  itself  is  not  a 
constant,  which  creates  classical  charge-density 
fluctuations  and  field  domains  in  MQWs.  For  non¬ 
steady  state,  the  initial  condition  for  Eq.  (3)  can  be 
set  as  Pj( 0)  =  e«2D  if  the  ac  electric  field  is  applied 
to  the  sample  after  t  =  0.  Eqs.  (3)-(5)  together 
(totally  2N  +  1  equations)  allow  us  to  simulta¬ 
neously  solve  for  the  charge-density  distributions 
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pk(t)  (or  chemical-potential  distributions  pk(t))  for 
k  =  1,2, . . .  ,N,  as  well  as  for  the  local-field  dis¬ 
tributions  Sk{t)  for  k  =  0, 1,2, . . .  ,7V  beyond  steady 
state.  For  steady  state  with  dpj{t)/dt  =  0,  only 
forward  contributions  will  stay.  In  this  case,  we 
only  need  to  replace  Eq.  (3)  at  each  moment  by  [8] 

JjD[Vp  Pj+U^j]  =  Mo  A] 

with  j  =  1, 2, . . .  ,N.  Under  this  condition,  the 
measured  tunneling  current  density  simply  equals 
J^[p0, pl,S’o\.  But  this  does  not  imply  a  uniform 
distribution  of  charge  densities  and  electric  fields. 


3.  Quantum  charge  fluctuations 

It  has  been  known  for  a  long  time  that  a  uni¬ 
form-field  distribution  will  underestimate  or 
overestimate  the  sequential-tunneling  current  in 
MQWs  at  high  or  low  voltages,  respectively,  [8] 
while  the  charge  density  will  remain  in  its  equi¬ 
librium  value  even  when  the  uniform  electric  field 
is  an  ac  field.  In  the  classical  field-domain  model 
[6,8],  quantum  properties  of  electrons  in  QWs  have 
been  ignored.  On  the  other  hand,  it  has  also  been 
shown  that  the  quantum  non-adiabatic  effects  give 
rise  to  a  residual  current  at  zero  of  applied  ac 
voltage  due  to  a  space-charge-field  effect  in  the 
presence  of  a  uniform  ac  electric  field  [2,5,11,12]. 
Therefore,  it  is  very  important  to  include,  simul¬ 
taneously,  both  the  classical  field-domain  effect 
and  the  quantum  non-adiabatic  effect  on  the  se¬ 
quential  tunneling  of  electrons  in  MQWs  when  an 
ac  electric  field  is  applied  to  the  sample. 

When  the  quantum  non-adiabatic  effect  is  in¬ 
cluded,  Eq.  (3)  for  j  =  1, 2, ...  ,7V  should  be  mod¬ 
ified  to 

dPj(t) 

d  t 

=  Pc(  1  -  PoY~lJ^\Po,  Ml  «>  *o(0M*o(0] 
+Ec^(i-Ecyi-1 

k=  1 

X  JkD\llk(t)  -  eLB^(t),pk+l(t),S”k(t)} 
x  e[^k(t)\-Pc(l-Pc)N-jJ^[p0,pN(t)^N(t)} 


x0[-£N(t)]-PcJ2(l-p°)k~J-1 

k=j+ 1 

x  0[fj( t)\  +jf[pj{t)-eL^\t),pj_M^_M 
Lu  d  £™{i) 

x  0[—Sj_x{t)\  — -^CQW  j  (6) 


where  the  last  term  represents  the  non-adiabatic 
increase  of  charge  density  in  the  y'th  QW,  is 
the  non-adiabatic  space-charge  field  in  the  yth  QW 
and  barrier,  is  the  cross-sectional  area  of  the 
sample  and  CQW  =  {m*e2£f /nJi2)f0{E0  —  pe(n2r>,  T)] 
is  the  quantum  capacitance  with  chemical  poten¬ 
tial  tk(w2D,  T)  for  an  equilibrium  two-dimensional 
electron  gas  in  QWs.  The  non-adiabatic  space- 
charge  fields  S™(t)  for  j  =  1 , 2, . . . ,  TV  in  Eq.  (6)  are 
determined  by  the  following  differential  equations 
[5,11,12] 


,  „  d Sj{t) 

qw^^-Cqw^t 


(7) 


which  contains  a  “quantum  displacement”  current 
due  to  Cqw  as  a  source  term.  The  non-adiabatic 
change  of  current  density  in  Eq.  (7)  is 

A J™(t) 

-  Jf  [/ Uj(t )  -  eLB<$™ ( t ) ,  A*y_,(0 ,  ,  {tm-sy  1  (/)] 

-  JT  [My(0 >  Pj+ 1  (0 ,  0[#j{t)] 

with  (t)  =  0.  The  charge-density  fluctuations 
in  Eq.  (6)  now  contain  both  adiabatic  and  non- 
adiabatic  contributions.  The  last  term  in  Eq.  (6) 
represents  the  non-adiabatic  contributions.  The 
total  field  S'k(t)  and  the  average  field  Sj{t)  in  Eqs. 
(6)  and  (7)  are  defined  by 


'm' 

.  \&k(t)  +  <nt-i(7)]/2_ 

Simultaneously,  Eq.  (4)  should  also  be  modified  to 

g^t)- g^t) 

=  —  -  en2D\, 

£oer 


(8) 
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where  C0  =  e0er^/ZB  is  the  classical  geometric 
capacitance.  For  non-steady  state,  the  total  mea¬ 
sured  current  density  for  Vb(t)  ^  0  is  given  by 

•C(+)M 

=  (1  -Pefjg{H,,MQ,*o(OM*o(t)\ 

+  52  ( 1  -  Po)N~kJlD[nk(t)  -  eL*gf(t) ,  pk+1  {t),S'k(t)] 

k=\ 

x  0[$'k(t)\  +y°[p0,pN(t),S:N(t)\6[-S'N(t)\ 
LBC0dSN(t) 

sr  d  t  ’  [  ] 

which  includes  both  conduction  and  dielectric 
displacement  currents,  where  the  last  term  repre¬ 
sents  the  displacement  current  from  the  geometric 
capacitance.  In  this  case,  the  in-flowing  injection 
current  is  not  equal  to  the  out-flowing  conduction 
current  from  the  sample  due  to  dpj(t)/dt  ^  0.  On 
the  other  hand,  for  Vb(t)  <  0  we  have 

■CHM 

+  E(i  -PcV 

k=  1 

X  42Dk(0  -  eUSf{i),pk_ i(t),4_i(t)] 
x  +  J^[pbl  p{{t),So{t)]d[S0{t)\ 

LjjCp  dtfpQ) 

y  d?  1  ’ 

in  steady  state,  however,  there  is  no  non-adiabatic 
charge-density  fluctuation,  and  the  total  measured 
current  density  is  mainly  determined  by  the  injec¬ 
tion  current,  which  is  given  for  Vb(t)  ^  0  by 

(id 


and  can  be  modified  by  the  local  field  at  the  emitter 
barrier.  Moreover,  due  to  dpj(t) /dt  ^  0  in  non¬ 
steady  state  we  can  define  a  total  differential  ca¬ 
pacitance  for  the  MQW  structure 


cdi(t)  =  jrCi(t) 


j= * 


dVbjt) 

dt 


h  *  ’ 

(12) 


which  is  time-dependent  and  different  from  both 
Co  and  Cqw- 


From  Eq.  (8)  we  know  that  both  the  field- 
domain  and  non-adiabatic  effects  will  cause 
charge-density  fluctuations  in  QWs  in  the  presence 
of  an  ac  electric  field.  The  quantum  capacitance 
Cqw  only  enters  into  Eq.  (6)  for  the  charge-density 
fluctuations  but  not  into  Eq.  (9)  for  the  total 
measured  current  density  J“(±)(t).  Instead,  the 
conduction  current  is  modified  by  the  non- 
adiabatic  space-charge  field  $y(t)  which  is  in¬ 
duced  by  the  “quantum  displacement”  current  as 
shown  in  Eq.  (7).  On  the  other  hand,  the  geometric 
capacitance  C0  directly  modifies  the  total  measured 
current  in  Eq.  (9)  as  a  contribution  from  the  di¬ 
electric  displacement  current  but  does  not  enter 
into  the  charge  fluctuations  in  Eq.  (6).  When  only 
the  forward  contributions  are  included,  the  sum  of 
dielectric  displacement  and  conduction  currents 
flowing  into  and  out  of  a  QW  is  a  constant  for 
Pc=  1.  Because  the  conduction  current  flowing 
into  the  first  QW  is  simply  the  injection  current 
from  the  top  contact  layer,  the  change  in  the 
conduction  currents  flowing  into  different  QWs  is 
determined  by  the  variation  of  the  dielectric  dis¬ 
placement  currents  due  to  the  non-uniform  elec¬ 
tric-field  distribution  inside  the  whole  system. 
From  Eq.  (7)  we  further  find  that  even  under  a 
uniform  ac  electric  field,  the  non-adiabatic  con¬ 
duction  current  density  flowing  through  each  QW 
is  not  equal  to  the  sum  of  adiabatic  sequential- 
tunneling  current  flowing  out  of  the  QW  and 
the  “quantum  displacement”  current  density 
L^Cowdig j(t) / dt  because  d<Sy{t)/dt  ^  0. 


4.  Non-adiabatic  self-consistent  Hartree  model 


The  non-equilibrium  electron  distribution 
function  in  a  shifted  Fermi-Dirac  model  [13]  can 
be  written  as 


/"»(k)=/'i0(%+Ak|  +  Ai7,), 
where  C(Ek)  is  defined  by 

cm  =  1 1  +  exp 


Ek  —  /(0(«2D,  T) 
hT 


(13) 


(14) 


For  the  shifted  Fermi-Dirac  model  in  Eq.  (13), 
there  exists  a  local  charge-density  fluctuation  for 
each  electron  state  |k),  defined  by 
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jSp*(t)  =  ^\roQO-fp(E*)] 

e/dAEk\\  8/HflOl 

r\  At  )X  m  J’  [  ’ 

where  Y  is  the  volume  of  the  sample. 

For  a  MQW  system  with  thick  barrier  layers, 
the  adiabatic  sequential-tunneling  current  density 
flowing  in  the  z  direction  (growth  direction  and 
perpendicular  to  the  QW  layers)  is  found  to  be 

Jn{t) 

=  %  E  I ^ |] [ r  (Ek)  -/"» ( Ek  +  e\fib |iB)] , 

7  k 

(16) 

where  v\  is  the  group  velocity  of  electrons  in  the  z 
direction,  and  8b  stands  for  the  applied  electric 
field.  If  C2zt  A>  1  with  t,  being  the  electron  se¬ 
quential-tunneling  time,  Y[Ek,  \8b\]  has  to  be 
found  by  solving  a  time-dependent  Schrodinger 
equation.  Otherwise,  Y[Ek,  |(fb|]  can  be  calculated 
from  a  static  Schrodinger  equation  at  each  time  t  if 
Qxt  <C  1.  We  will  be  only  interested  in  the  latter 
case  with  Qrt  <C  1  hereafter.  If  we  replace  fll"(Ek) 
to  leading-order  approximation  by  the  equilibrium 
value  fg°(Ek)  in  Eq.  (14)  for  faster  electron  energy 
relaxation  processes  due  to  inelastic  scattering  of 
electrons  compared  to  the  electron  sequential 
tunneling,  and  replace  the  electron  group  velocity 
tfk  by  a  drift  velocity  ud[<?b]  (a  statistically-averaged 
group  velocity)  of  electrons  in  a  bulk  material,  Eq. 
(16)  reduces  to  Levine’s  sequential-tunneling 
model  [1] 

b]  =  E  f[Ek,  Kbl  }[fo^(Ek) 

'  k 


—  fo°(Ek  +  e\Sb\L^)\,  (17) 

where  vd[Sb]  =  ( ezv/m*)Sb ,  and  the  momentum- 
relaxation  time  tp  is  given  by 


In  Eq.  (17),  J^[Sb]/evA[Sb]  can  be  equivalently 
viewed  as  a  three-dimensional  tunneling-electron 
density  which  depends  on  Sb,  T  and  h2d- 

From  now  on,  we  limit  ourselves  to  an  electri¬ 
cal-quantum  limit  where  only  the  ground  subband 


of  the  quantum  well  is  occupied  by  electrons  at  low 
temperatures  and  low  electron  densities.  The 
electron  kinetic  energy  of  the  ground  subband 
(measured  from  the  edge  E0)  is  given  by  Ek  = 
h2k2 /2m* .  In  the  current-surge  model  [5,11,12],  we 
assume  that  A Ek  is  associated  with  the  global 
fluctuation  (independent  of  individual  electron 
state)  of  the  chemical  potential  of  electrons  in  the 
QW.  By  writing  AEk  =  —  A/x  =  fi0(n2D,  T)  —  n{t) 
for  the  global  chemical-potential  fluctuation, 
where  / i(t )  and  n0(n2r>,T)  are,  respectively,  the 
transient  chemical  potential  for  electron  density 
ne(t)  and  that  for  an  equilibrium  electron  gas  in 
quantum  wells,  we  get 

dAEk  _  dAEk  0/i  d«e 

d t  ~  dt  0«e  dt  '  1  ’ 

We  further  introduce  a  spatially-averaged  non- 
adiabatic  space-charge  field  Sm(t)  which  is  defined 
by  [5,11,12] 

A Ek  =  e£m(t)LB,  (20) 

where  £’m(t)  measures  the  non-adiabatic  reduction 
of  the  electron  chemical  potential  in  QWs  under  a 
uniform  applied  electric  field.  If  we  use  Levine’s 
sequential-tunneling  model  in  Eq.  ( 1 7),  we  find  the 
following  non-adiabatic  current  density  due  to  the 
existence  of  this  non-adiabatic  space-charge  field 

«?na  (t) 

A Jna(f)  -  J*>-**'-  \Sb  +  <?na]  -  J ^  [rf’b] ,  (21) 

where  Jm' [S'b]  has  been  given  in  Eq.  (17).  In  Eq. 
(21),  the  first  term  can  be  viewed  as  an  equivalent 
capture  current  into  the  QW,  while  the  second 
term  can  be  regarded  as  a  sequential-tunneling 
current  flowing  out  of  the  QW. 

For  a  QW,  the  electron  density  will  be  constant 
if  the  conduction  currents  flowing  in  and  out  of  the 
well  are  equal.  The  variation  of  the  charge  density 
in  the  well  is  created  by  an  imbalance  in  conduc¬ 
tion  currents.  The  charge-current  conservation  law 
requires 

rjhp{t)  EE  Yjt  E 5 PM  =  YAJUt)-  (22) 

The  left-hand  side  of  Eq.  (22)  represents  the  non- 
adiabatic  charge  increase  inside  the  well,  while  the 
right-hand-side  of  the  equation  stands  for  the  net 
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increase  of  charge  due  to  non-adiabatic  current 
flowing  into  the  QW. 

If  Qrt  <C  1,  the  ground-state  electron  wave 
function  cj)]  (z,  t)  inside  the  QW  within  the  self- 
consistent  Hartree  model  is  determined  by  [14] 


fr  d  / 

'  1  d  \ 

2  dz  \ 

v  m*(z)  dz  J 

-  eSb(t)z  +  UqW{z ) 


+  + h(z,  t) 


4>i(z,  t)  =  E0(t)<l>l(z,t), 


(23) 


where  E0(t)  is  the  time-dependent  ground-subband 
edge,  the  electron  effective  mass  m*  (z)  is  ww  in  the 
well  and  mb  in  the  barrier,  and  t/QW(z)  is  zero  in¬ 
side  the  well  but  V0  outside  the  well.  Within  the 
adiabatic  limit,  we  have  p(t )  =  p0(n2D,T),  other¬ 
wise  S p(t)  ^  0  for  non-adiabatic  cases.  The  Har¬ 
tree  potential  tT'H(z,  t)  in  Eq.  (23)  can  be  found 
from  the  Poisson  equation 


d_ 

dz 


Cr(Z)  H(Z’  ^ 


=  —[N+z)  -«e(z,  0]. 
^0 


(24) 


where  donors  are  assumed  completely  ionized,  and 
the  relative  dielectric  constant  er(z)  is  ew  in  the  well 
and  eB  in  the  barrier.  N+z)  in  Eq.  (24)  is  the  static 
profile  of  donor  doping  for  the  single  quantum 
well,  ne(z,t)=  1 (z, 0|2«e(0  is  the  density  func¬ 
tion,  and 


r+oo 

«e(t)  =  «2D  +  § P(t)  =  «2D  +  P2D  /  dE6f(E,  t)  , 

Jo 

(25) 

where  p2D  =  (m^/nfr)  is  the  density  of  states  for 
two-dimensional  electrons  in  the  QW,  5 f(E,t) 
represents  the  local  non-adiabatic  fluctuation  of 
the  electron  distribution  function  in  energy  space. 
Here,  the  number  of  electrons  in  the  quantum  well 
is  not  a  constant  due  to  the  non-adiabatic  current 
flowing.  Moreover,  we  find  from  Eqs.  (15),  (19) 
and  (25) 

jpp(') 


=  ep2  D 


dE-8f(E,t) 


d  r 

e2Lbp2D—£’b(t)  J 


r+oo 

/  d  E 

(E)  1 

lo 

dE 

(26) 


Applying  Eq.  (22)  and  using  Eqs.  (21)  and  (26),  we 
find  the  following  integral  equation  for  8 f(E,  t )  by 
using  Levine’s  model  in  Eq.  (17) 


eP  2D 


dE-hf{E,t) 


-e  Lbp2D 
ePlD 


d<$b{t) 

dr 


r+oo 

/  dE 

dfo°  (E)  1 

Jo 

dE 

Ew 

r+oo 


{  ^  b]  +  St’d  [8/] } 


x  dE^r[E+E0M,rH}x[f^(E) 


+  8f(E,t)-f«>(E+e\gb\LB)-bf(E+e\gb\LB,t)\ 

r+oo 


ePlD 


Vd  [gb 


dE:T[E+E^\\gb\-r^}[C{E) 


(0)i 


-ff{E+e\Sb\L*)]=0, 


(27) 


where  +\E  +  Eq,  \gb\\ f''H]  is  the  quantum  trans¬ 
mission  of  electrons  and  can  be  calculated  by  it¬ 
eration  (see  Appendix  A).  In  Eq.  (27),  i'+(z,  t)  and 
E0(t)  are  written  simply  as  tE'h  and  E0.  The  adia¬ 
batic  quantities  +h(z,  t)  and  E(0°\t )  can  be  ob¬ 
tained  by  simply  setting  8 p(t)  =  0  in  Eq.  (24)  and 
+h(z,  t )  =  rf{z,  t)  in  Eq.  (23).  Moreover,  8t>d[§/] 
introduced  in  Eq.  (27)  is  calculated  to  be 


§M5/]  =  - 


(28) 


Finally,  Eq.  (27)  leads  us  to  the  non-adiabatic 
dynamical  differential  equation  for  §/(£,?) 


^8 f{E,t)-eLb 


d^b(Q  d/o*0 (E) 

dr 


S E 


-4-{fdKb]  +  8rd[5/]}^[E+Eo,Kb|;^H]  x  [ ff(E) 

Ly/ 

+  5/(£,r)  —  fo°  (E  +  e\g  b\Lb)  —  d>f(E  +  e\Sb\LBp)\ 


+  y— t’a  [gb]3r[E + 40) .  Kbl ;  r^]  [/*»  (e) 

-f«>(E  +  e\gb\LB)]=0, 


-(0)1 


(29) 


where  the  initial  condition  is  chosen  to  be 
8/(£,  t)  =  0  at  t  =  0  if  the  ac  electric  field  is  ap¬ 
plied  to  the  system  after  t  =  0.  8/(£,  t )  has  a  lower 
bound  which  is  determined  by  the  condition 

8f(E,t)+fS°(E)=0. 

For  small  A p,  the  first  term  in  Eq.  (29)  can  be 
approximated  to  leading  order  by 
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0  „,r  ,  0A/( 

—  8f(E,  t)  «  —A 
dt  y  J  dt 


g/T(g) 

dE 


(30) 


Similarly,  the  third  term  in  Eq.  (29)  can  be  ap¬ 
proximated  as 

*■[£  +  Eo,  Kbl;  ru\lC(E)  +  8f(E , 0 
-/0"°(^  +  e|^b|iB)  -  6/(£+  e|<?b|W)] 

»  +  4°\  I  Al;  <0)  +  (§<)/5«2d)p2DA/i] 

x  [C+A"(£)  ™/0"“+a"(£  +  e|*b|IB)],  (31) 

where  (SiG^/S^d)  =  (e2/2e0ew<7TF)  and  </TF  = 
(e2/2e0ew)p2D  from  the  Thomas-Fermi  model  [14]. 
By  recalling  A/t  =  -eSna(t)L^,  Eq.  (29)  results  in 
the  current-surge  model  [11],  where 


Cqw  =  p2D 

Jo 


p  +oo 

/  d  E 

0/;°(^)i 

Jo 

dE 

J^[S  b 


ePlD 


Vd  [<*  b 


d  EST 


JW  JO 

_  /■«> 


£  +  £f,|<fb|;^ 


(0) 


X  [/H^-Z^  +  W*)], 


JT 


£  +  £ 


(0)  | 
0  5  I 


fb|;  ~  e^naAfi 


E  +  E{f\Sb  +  Sm\-r^ 


(0) 


and  vd[Sb]  +  8rd[8/]  «  vA[Sh  +  <fna]  if  we  set  tp  « 
Ln/vy  with  vF  being  the  electron  group  velocity  at 
the  Fermi  energy.  The  non-adiabatic  space-charge 
field  (?na (f)  introduced  by  Ap  =  -e<f>m(t)Lv  can  be 
calculated  from 


i  p+oo 

<?na  (t)  =  ~—  dESf(E,t),  (32) 

eLb  J o 

which  becomes  positive  if  8 p(t)  <  0. 


5.  Numerical  results 

In  this  section,  we  first  present  numerical  results 
for  distributions  of  both  local  fields  and  charge 
densities  in  the  presence  of  an  applied  ac  electric 
field,  including  effects  of  classical  and  quantum 
charge  fluctuations.  After  this,  we  present  results 
for  fluctuating  electron  distributions  and  Hartree 
potentials.  The  voltage  is  defined  to  be  Vb (t)  = 


&b{f)Lt  with  $b(t)  =  <odc  +  <?ac(t)  and  Sac(t)  = 
S'm  s\n{2nt/Tp)  for  t  ^  0. 

The  sample  1  we  consider  for  showing  domain 
effects  is  an  AlGaAs/GaAs  MQW  structure.  The 
total  number  of  QWs  is  N  =  10,  with  eleven  bar¬ 
riers.  The  parameters  for  sample  1  are:  well 
thickness  Z.w  =  75  A,  barrier  thickness  Lb  =  339  A, 
barrier  height  V0  =  224.5  meV,  electron  effective 
mass  m*  =  0.065»?0  with  free  electron  mass  mb, 
electron  areal  density  «2d  =  5  x  10u  cm-2,  con¬ 
tact-layer  electron  concentration  nc  =  6  x  1017 
cm-3,  cross-sectional  area  =  10  4  cnr,  capture 
probability  Pc  =  0.5,  saturation  velocity  vs  = 
2  x  106  cm/s,  saturation  field  <Ss  =  2  kV/cm,  and 
relative  dielectric  constant  er  =  12.  For  this  sam¬ 
ple,  the  ground-state  subband  edge  is  calculated  to 
be  E0  =  44.1  meV.  For  the  applied  ac  electric  field, 
S’m  =  5  kV/cm,  <fdc  =  0  and  Tp  =  0.1  s.  The  sample 
2  we  chose  for  the  non-adiabatic  self-consistent 
Hartree  model  is  also  an  AlGaAs/GaAs  MQW 
structure.  The  parameters  for  sample  2  are: 
Lw  =  80  A,  Lb  =  300  A,  bb  =  331  meV,  electron 
effective  mass  in  well  mw  =  0.067m0,  electron  ef¬ 
fective  mass  in  barrier  mw  =  0.092»?0,  «2d  = 
4  x  1011  cm-2,  ST  =  2.25  x  10  4  cm2,  ds  =  2x  106 
cm/s,  $$  =  2  kV/cm,  dielectric  constant  in  well 
ew  =  12.0,  and  dielectric  constant  in  barrier 
ew  =  11.2.  Eq0)  is  calculated  to  be  44.5  meV.  For 
the  applied  electric  field,  <fdc  =  0.05  kV/cm  and 
Tp  =  4  s. 

We  show  in  Fig.  2  the  local  fields  S'j(t)  -  Sb(t) 
in  different  layers  of  sample  1  (in  (a))  and  the 
density  fluctuation  nj(t)  —  m2d  in  different  QW’s  (in 
(b))  at  several  times  t/Tp  for  T  =  77  K.  From  (a) 
we  find  that  the  field-domain  effect  is  negligible 
at  f/7p  =  0.05  (very  small  applied  field)  because 
both  the  injection  current  and  the  sequential- 
tunneling  current  JA2D  are  both  extremely  small  in 
this  case.  With  the  increase  of  Sb{t),  i.e.  t/Tp  in¬ 
creases  from  0.05  to  0.25,  fields  close  to  the  emitter 
barrier  are  enhanced  dramatically  relative  to  the 
uniform  field  Sb  (t).  This  is  a  result  of  the  huge 
current  imbalance  J;345  -C  Jfn  under  the  uniform 
field  Sb{t),  as  explained  in  Fig.  1(b).  At  the  same 
time,  fields  close  to  the  receiver  barrier  are  sup¬ 
pressed  almost  to  zero.  From  (b)  we  find  that  when 
Sb(t)  is  large,  densities  close  to  the  emitter  barrier 
are  greatly  reduced  with  respect  to  the  equilibrium 
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Fig.  2.  Calculated  local  fields  Sft)  —  $b{t)  in  (a)  in  different  barrier  layers  and  density  fluctuations  nft)  —  n2 d  in  (b)  inside  different 
QWs  of  sample  1  at  times  t/Tv  =  0.05,  0.11,  0.20  and  0.25  s.  The  parameters  used  in  calculations  are  given  in  the  text. 


value  «2D.  This  is  accompanied  by  a  great  en¬ 
hancement  of  the  local  field  SQ(t)  »  <Sh{t)  at  the 
emitter  barrier,  as  shown  in  (a).  However,  densities 
close  to  the  receiver  barrier  remain  near  m2D  due  to 
the  suppressed  local  field  S‘N(t)  «  0.  It  is  obvious 
from  (b)  that  some  electrons  have  been  removed 
from  the  sample  since  V  ■  [«/(*)  —  »2d]  <  0,  which 
is  true  even  for  steady  state.  The  calculation  done 
here  corresponds  to  a  non-steady  state.  Therefore, 
the  net  number  of  electrons  removed  from  the 
sample  changes  with  time  under  an  ac  voltage, 
leading  to  a  differential  capacitance  CJdi(t)  (see  Eq. 
(12)). 


Fig.  3  compares  local  fields  Sj(t)  —  Sb(t)  in 
sample  1  at  t/Tp  =  0.25  as  a  function  of  barrier 
index  j  for  different  values  of  <?m,  Tp,  T  in  (a)  and 
different  values  of  N,  n2 d,  nc  in  (b).  From  (a)  we 
find  that  the  field-domain  effect  is  negligible  at 
7X65  K  due  to  very  small  injection  and  sequen¬ 
tial-tunneling  currents  at  these  temperatures.  The 
bigger  the  field  amplitude  Sm  is,  the  larger  S0(t) 
will  be.  A  smaller  Tp  leads  to  a  negative  SN{t)  on 
the  receiver  barrier  due  to  the  strong  non-steady 
effect.  This  is  completely  different  from  the  steady- 
state  results  [8]  in  which  S'N{t)  will  always  be 
positive.  Furthermore,  we  find  from  (b)  that  the 


Fig.  3.  Calculated  local  fields  $  ft)  —  (t)  in  different  barrier  layers  of  sample  1  at  the  time  t/Tp  =  0.25  s.  Here,  the  comparison  of 
$  ft)  —  Sh (t)  with  different  values  of  Tp  and  T  are  presented  in  (a)  and  the  comparison  of  those  with  different  values  of  N,  n2 d  and 

nc  are  shown  in  (b).  The  changed  parameters  are  indicated  in  figures.  The  other  parameters  used  in  calculations  are  the  same  as  those  in 
Fig.  2. 


D.  Huang,  D.A.  Cardimona  /  Infrared  Physics  &  Technology  44  (2003)  487-501 


497 


5.035 


5.030 


5.025 


5.020 

0  2  4  6  8  10 

Well  Index  j 

Fig.  4.  Calculated  non-adiabatic  space-charge  fields  <^"a( t )  in 
sample  1  at  t/Tp  =  0.25  as  a  function  of  well  index  j  for  T  —  40 
and  77  K  (left  scaled),  respectively.  The  other  parameters  used 
in  calculations  are  the  same  as  those  in  Fig.  2. 
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smaller  the  number  of  QW’s  N  is,  the  lower  <f0(f) 
is.  The  increase  of  n2 d  causes  a  larger  sequential- 
tunneling  current,  leading  to  a  larger  value  of  S'o(t) 
due  to  an  enhanced  current  imbalance  between  J® 
and  JfD.  Conversely,  the  increase  of  nc  introduces 
a  bigger  injection  current,  leading  to  a  smaller 
value  of  S\)(t)  due  to  a  suppressed  current  imbal¬ 
ance  between  and  J/D. 

Effects  of  classical  charge  fluctuations  depend 
on  the  geometric  capacitance  C0,  as  shown  by  Eq. 


(4).  However,  the  quantum  capacitance  CQW  starts 
to  play  a  role  when  the  non-adiabatic  effect  is  in¬ 
cluded,  as  seen  from  Eq.  (8).  The  striking  thing  is 
that  the  effects  of  quantum  charge  fluctuations 
become  more  and  more  important  as  T  C  65  K, 
while  the  effects  of  classical  charge  fluctuations  are 
negligible  at  these  temperatures.  We  compare  in 
Fig.  4  two  calculated  non-adiabatic  fields  <?“(/)  at 
t/Tp  =  0.25  in  sample  1  as  a  function  of  well  index 
j  for  T  =  40  and  77  K.  From  the  figure  we  know 
that  (f“a(f)  decreases  with  j  due  to  field-domain 
effects  at  T  =  77  K.  However,  <f”a(t)  becomes  in¬ 
dependent  of  j  at  T  =  40  K.  More  importantly, 
<f”a(t)  increases  with  reducing  T. 

The  effect  of  non-adiabatic  field  <?"a(t)  can  be 
seen  more  clearly  from  the  calculated  local  fields 
S'j(t)  —  <fb(f)  in  different  barrier  layers  of  sample  1 
in  Fig.  5(a)  and  density  fluctuations  rij(t)  —  w2d  in 
different  QW’s  in  Fig.  5(b)  at  t/Tp  =  0.25  with 
various  values  of  S'm.  7j,  and  7’.  By  comparing  solid 
and  dashed  curves  in  (a)  we  find  that  S0(t)  is  re¬ 
duced  by  a  factor  of  2  as  Sm  =  5  kV/cm,  Tp  =  0.1  s 
and  T  =  77  K  (squares).  Moreover,  S0(t)  decreases 
even  more  due  to  (f”a(f)  when  7j,  is  reduced  to  0.05 
s  (stars)  due  to  stronger  non-steady  effects,  but  it 
decreases  much  less  when  T  is  reduced  to  40  K 
(triangles)  due  to  smaller  injection  and  sequential- 
tunneling  currents.  However,  S0(t)  is  enhanced 
by  non-adiabatic  effects  when  Sm  is  reduced  to 


Fig.  5.  Calculated  local  fields  Sj(t)  —  $b(t)  in  different  barrier  layers  of  sample  1  and  density  fluctuations  nj{t)  —  «2d  inside  different 
QWs  at  t/Tp  =  0.25  with  different  values  of  (squares  and  circles),  Tp  (squares  and  stars)  and  T  (squares  and  triangles).  Sj{t)  —  <fb (0 
are  compared  in  (a)  with  (solid  curves)  and  without  (dashed  curves)  i f”a(t),  and  the  similar  comparisons  of  nj(t)  —  n2 d  are  presented  in 
(b).  The  changed  parameters  are  indicated  in  figures.  The  other  parameters  used  in  calculations  are  the  same  as  those  in  Fig.  2. 
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1  kV/cm  (circles).  On  the  other  hand,  from  (a)  we 
also  find  that  <SN{t)  becomes  a  much  more  negative 
value  for  Sm  =  5  kV/cm,  Tv  =  0.1  s  and  T  =77  K 
when  is  included  in  our  calculations  (solid 

curve  with  squares)  than  when  it  is  excluded  (da¬ 
shed  curve  with  squares).  The  same  situation  oc¬ 
curs  when  either  Tv  is  reduced  to  0.05  s  (stars)  or 
dfm  is  reduced  to  1  kV/cm  (circles)  but  with  a 
smaller  overall  magnitude  compared  to  the  curves 
with  the  squares.  The  features  observed  for  SN(t) 
in  Fig.  5(a)  will  affect  the  conduction  current  de¬ 
tected  at  the  receiver  layer.  Effects  of  quantum 
charge  fluctuations  are  reflected  in  the  calculated 
Hj(t)  —  «2d  in  (b),  where  the  reduction  of  charge  in 
the  QWs  is  greatly  increased  except  for  the  case 
with  small  value  of  <?m  (circles).  Although  there  is 
a  strong  dependence  of  Sj{t)  on  index  j  near  the 
emitter  barrier,  there  is  very  little  dependence  of 
rij(t)  on  j  there  after  $™(t)  is  included  in  the  cal¬ 
culations. 

In  order  to  gain  further  information  about  the 
local  change  in  the  non-adiabatic  electron  distri¬ 
bution  function,  we  display  5/(7?,  t)  of  sample  2  in 
Fig.  6  at  t/Tp  =  0.25  with  uniform  doping  for 
different  values  of  T.  From  the  figure  it  is  clear  that 
5/(7?,  t)  always  shows  a  negative  minimum  at 
llo(n2D,T)  since  it  is  proportional  to  cfp'{E) /dE. 
Since  the  Fermi  surface  broadens  with  increasing 
T,  we  find  from  the  figure  that  the  negative  mini¬ 
mum  is  partially  suppressed  and  broadened  (solid 
curve)  when  T  =  40  K  compared  to  that  at 


Fig.  6.  Calculated  non-adiabatic  change  in  distribution  func¬ 
tions  6 f(E,t)  at  t/Tp  =  0.25  for  electrons  in  sample  2  with 
uniformly  doped  QWs.  In  this  figure,  we  set  <?m  =  1  kV/cm  with 
T  =  40  K  (solid  curve)  and  T  =  20  K  (dashed  curve). 

T  =  20  K.  We  also  find  that  the  negative  minimum 
is  enhanced  when  Sm  is  increased  (not  shown). 

Fig.  7(a)  and  (b)  present  non-adiabatic  charging 
effects  in  the  uniformly-doped  QW  on  the  Hartree 
potentials  for  sample  2  at  T  =  40  K  and  Sm  =  1 
kV/cm.  From  (a)  we  find  that  the  positive  peak  in 
the  adiabatic  Hartree  potential  V^°\z,t)  at  the 
center  of  the  QW  is  greatly  suppressed  by  the  non- 
adiabatic  effect  at  t/Tp  =  0.25,  leaving  two  positive 
spikes  at  the  edges  of  the  QW.  Fig.  7(b)  shows  the 
comparison  between  non-adiabatic  Hartree  po- 


Fig.  7.  Calculated  position  z  dependence  of  non-adiabatic  (solid  curve)  and  adiabatic  (dashed  curve)  Flartree  potentials  of  sample  2  in 
(a)  at  T  =  40  K  and  =  1  kV/cm  with  uniform  doping  inside  the  QW  and  non-adiabatic  Hartree  potentials  FH(z,  t)  in  (b)  for 
t/Tp  =  0.25  (solid  curve)  and  t/Tp  =  0.75  (dashed  curve). 
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Fig.  8.  Calculated  non-adiabatic  electron  distribution  functions  /(£,  t)  in  (a)  and  logarithm  of  absolute  value  of  total  non-adiabatic 
sequential-tunneling  current  /na(f)  =  [Sb]  +  A/na(t)]  of  sample  2  as  a  function  of  applied  ac  electric  field  in  (b)  for  T  =  40  K 
and  Sm  =  1  kV/cm.  In  (a),  we  plot/(E,  t)  at  t/Tp  =  0.25  (dotted  curve)  and  t/Tp  =  0.75  (dashed  curve).  The  time-independent  adiabatic 
electron  distribution  function  (solid  curve)  is  also  shown  for  the  comparison.  In  (b),  the  currents  calculated  from  adiabatic 

(lower  curve)  and  non-adiabatic  (upper  curve)  electron  sequential  tunneling  are  compared  with  each  other. 


tentials  when  electrons  in  the  QW  are  either  re¬ 
moved  (t/Tp  =  0.25,  solid  curve)  or  added  (t/Tp  = 
0.75,  dashed  curve).  We  find  from  (b)  that  the  two 
positive  spikes  at  the  edges  of  the  QW  are  sup¬ 
pressed  but  two  negative  spikes  are  enhanced  when 
electrons  are  added  to  the  well. 

Finally,  we  display  in  Fig.  8  f(E,t)  of  sample  2 
at  t/Tv  =  0.25  (dotted  curve)  and  0.75  (dashed 
curve),  as  well  as  in  the  adiabatic  approximation, 
( f£°(E ),  solid  curve)  in  (a)  and  the  log10  |/na(0l  as  a 
function  of  Sb(t)  in  (b).  From  (a)  we  see  f(E,  t)  has 
a  fluctuation  in  time  with  respect  to  the  adiabatic 
fo°(E)  around  the  Fermi  energy.  Compared  with 
the  adiabatic  electron  sequential-tunneling  current 
(lower  solid  curve  with  A/na(f)  =  0)  in  (b),  the 
symmetry  of  log10  |/na(t)|  with  respect  to  the  posi¬ 
tive  (electrons  being  removed)  and  negative  (elec¬ 
trons  being  added)  peaks  of  Sb(t)  is  broken  in  the 
non-adiabatic  case  (upper  solid  curve).  A  small 
offset  [5]  of  the  non-adiabatic  log10  |/na(f)|  relative 
to  <Sb(t)  =  0  can  be  seen  by  comparing  upper  and 
lower  solid  curves. 


6.  Conclusions 

In  conclusion,  by  including  non-adiabatic 
space-charge-field  effects  we  have  generalized  the 
previous  theories  for  studying  field-domain  effects 


in  MQW  photodetectors  in  the  presence  of  an  ac 
voltage.  We  have  found  from  our  numerical  cal¬ 
culations  that  field-domain  effects  are  only  im¬ 
portant  at  high  temperatures  or  high  voltages, 
which  implies  the  existence  of  significant  injection 
and  sequential-tunneling  currents  in  the  system. 
We  have  further  found  that  non-adiabatic  effects 
become  much  more  visible  at  low  temperatures 
and  low  voltages  when  the  field-domain  effects  are 
negligible.  Furthermore,  we  have  derived  a  dy¬ 
namical  differential  equation  for  the  non-adiabatic 
electron  distribution  function  for  sequential- 
tunneling  current  flowing  through  a  MQW  system. 
Using  this  equation,  we  generalized  the  self-con¬ 
sistent  Hartree  model  for  the  calculation  of  non- 
adiabatic  electronic  states  in  a  quantum  well. 
Finally,  we  have  connected  the  current  quantum- 
statistical  theory  to  the  previously-proposed 
current-surge  model  with  a  leading-order  approx¬ 
imation. 

In  this  paper,  we  have  assumed  the  capture 
probability  is  independent  of  electric  field.  This  can 
be  justified  by  the  fact  that  the  capture  probability 
is  nearly  constant  at  low  electric  fields.  In  the 
presence  of  incident  photons,  the  conduction  cur¬ 
rent  flowing  through  the  MQW  sample  will  be  the 
sum  of  sequential-tunneling  and  photoexcited 
currents.  From  our  studies  in  this  paper,  we  predict 
that  the  field-domain  effects  which  are  significant  at 
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high  temperatures  or  high  photon  fluxes  will 
strongly  affect  both  tunneling-  and  photocurrents. 
On  the  other  hand,  non-adiabatic  effects  modify 
these  currents  at  low  temperatures  and  low  photon 
fluxes  when  the  current  flowing  through  the  MQWs 
is  small.  These  latter  conditions  are  of  utmost  im¬ 
portance  in  any  space-based  detector  applications. 
Furthermore,  only  the  self-consistent  Flartree 
model  is  employed.  The  exchange  interaction  be¬ 
tween  electrons  is  expected  to  be  very  small  [16]  at 
T  =  40  K  and  has  been  neglected. 

The  time  scale  for  observing  the  non-adiabatic 
space-charge  effect  requires 

t  <  RtCQw£f  <  tt  <C  271  /Q, 

where  Rt  =  {LB/ £R)[QJ,lv{Sb}/QS\,Yl  is  the  differ¬ 
ential  tunneling  resistance,  depending  on  and  T. 
Further,  Qrt  -C  1  ensures  the  electrons  see  only  an 
instantaneous  ac  electric  field  during  the  sequen- 


Appendix  A 


The  quantum-mechanical  transmission  coeffi¬ 
cient  3T[E  +  Eu$b\i'~ H]  used  to  evaluate  the  tun¬ 
neling  current  in  Eq.  (27)  can  be  found  from  the 
following  backward  iteration  [15]  at  each  time  t 


-i(0  —  {  2  + 


2  +  ^\uf-eSb{t){j-\)A 

Pa  L 


+  Vy{t)-E-El{t)\  \4>j{t)-4>j+l{t) 

(A.1) 


for  where  r/j/t)  =  r/qfz,-,?),  V"(t)  = 

VH (zj,  t),  Ed  =  h2/2mBA2,  A  =  LB/NB  and  NB  is  the 
number  of  slabs  (thickness  A)  within  the  barrier 
layer.  Flere,  Uf  =  0  for  j  =  0  and  j  =  NB  +  1 . 
Otherwise,  Uj  =  V0.  The  ending  boundary  condi¬ 
tion  of  Eq.  (A.l)  produces 


{1  -i  V 


1EA\E+E,(t)  -  U 


exp(i/VB&') 


(A.2) 


tial-tunneling  process.  Finally,  t  <  Rt  Cq\Vs/  en¬ 
sures  the  observation  of  a  non-adiabatic  space- 
charge  effect  inside  the  QW.  The  tunneling  time  tt 
can  be  estimated  from  rt  ~  e/[J,l°[S\,]£R}.  For  a 
superlattice,  we  take  Jli0[$h\sd  =  1  pA,  leading  to 
t,  =  0.1  ps  and  Q  <C  1013  FIz  from  Qrt  <C  1.  For  a 
MQW  system,  we  take  J^[£h\£R  =10  pA,  leading 
to  t,  =  10  ns  and  Q  <C  10s  Hz.  Difficulties  in  ob¬ 
serving  the  non-adiabatic  effect  may  come  from 
the  small  QW  capacitance  CQW^  ~  100  pF  in  the 
requirement  t  <  RtC()\V<f’.  For  a  superlattice,  we 
take  Rt  =  104  Q,  and  then  t  <  1CF6  s  is  required 
(very  hard  to  observe).  For  a  MQW  system,  on  the 
other  hand,  we  take  Rt  =  1011  Q,  which  implies 
t  <  10  s  (very  easy  to  observe). 


where  V  =  (A/h)^/2mB(E  +  e£b(t)LB).  From  the 
solution  of  Eq.  (A.l)  we  find  the  quantum- 
mechanical  transmission  of  electrons  from 


^[E  +  EuS^rn} 


1  jE  +  eS’b(t)LB 

jdf  V  E 


(A.3) 


where  |S|2  =  [|a|2  +  \b\2  +  2Re(ah*)]/4.  Here,  the 
two  complex  numbers,  a  and  b,  are  defined  by  the 
starting  boundary  condition  of  Eq.  (A.l) 


a 

b 

.-(i /2k){<t>2(t)  -  <t>o(t)) 

with  k  =  (A/h)  \J2mBE. 


(A.4) 
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